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I introduce a critical-state theory incorporating both flux cutting and flux transport to calculate 
the magnetic-field and current-density distributions inside a type-II superconducting cylinder at its 
critical current in a longitudinal applied magnetic field. The theory is an extension of the elliptic 
critical-state model introduced by Romero-Salazar and Perez- Rodriguez. The vortex dynamics 
depend in detail upon two nonlinear effective resistivities for flux cutting [p\\) and flux flow 
and their ratio r = P\\/p±- When r < 1, the low relative efficiency of flux cutting in reducing 
the magnitude of the internal magnetic-flux density leads to a paramagnetic longitudinal magnetic 
moment. As a model for understanding the experimentally observed interrelationship between the 
critical currents for flux cutting and depinning, I calculate the forces on a helical vortex arc stretched 
between two pinning centers when the vortex is subjected to a current density of arbitrary angle 
0. Simultaneous initiation of flux cutting and flux transport occurs at the critical current density 
Jc(0) that makes the vortex arc unstable. 

PACS numbers: 74.25.Sv,74.25.0p,74.25.Ha,74.25.Wx 



I. INTRODUCTION 

The behavior of type-II superconductors carrying cur- 
rent in a perpendicular applied magnetic field is well un- 
derstood in terms of the critical-state theory, first in- 
troduced by Bean.^'^ The fundamental idea is that the 
penetration of magnetic fields in the form of quantized 
vortices^ into practical type-II superconductors is im- 
peded by an array of pinning centers.^ The magnetic- 
fiux distribution can be put into a variety of metastable 
states, and changes in the fiux distribution can occur 
only when the magnitude of the local current density 
J = \/ X H exceeds the critical-current density Jc, the 
threshold for depinning and flux transport. The driving 
force per unit volume on a vortex array carrying mag- 
netic flux density B is F = J xB^ and the vortices move 
whenever the magnitude of this force locally exceeds the 
average pinning force per unit volume, Fp = JcB. Mo- 
tion of the vortices with a local velocity v gives rise to a 
local flux-transport electric field^ E = B x v perpendic- 
ular to B. This critical- state theory has served us well 
in providing a basis for practical applications of super- 
conductivity, such as in magnet technology,^ where the 
magnetic fields generated are practically all perpendicu- 
lar to the currents that generate the fields. In addition, 
this theory has permitted a good understanding of the 
ac losses in many electric power applications,^ since the 
local rate of power dissipation J • E is easily calculable 
when J = I J| is just above Jc- 

On the other hand, the corresponding theory for the 
behavior of type-II superconductors carrying current in 
a parallel applied magnetic field or in a field at an arbi- 
trary angle relative to the current is not well developed. 
With respect to the standard critical- state theory, two 
key questions are: (a) If J is not perpendicular to B^ 
how large can the (force-free) component of J parallel to 



B be? Is there any limit to this component below the 
depairing current density? (b) Superconducting wires 
subjected to a parallel magnetic field experimentally ex- 
hibit a finite electric field E with a component parallel to 
B when the current carried by the wire exceeds the crit- 
ical current.^ How is this component of the electric field 
produced? Flux cutting^^^ provides natural answers to 
both of these questions: When Jy (the component of J 
parallel to B) is small, E\\ (the component of J parallel 
to -B) is zero, but when the magnitude of Jy exceeds the 
threshold for flux cutting, flux-cutting processes initiated 
by local helical instabilities^"-^^ generate a finite value of 
E\\. 

Analogous processes occur in rotating superfluid ^He, 
where thermal counterfiow parallel to the vortices pro- 
duces turbulence initiated by the Glaberson-Donnelly 
helical instability.-^'' The resulting energy input to the 
vortex system is dissipated at the microscopic level by 
vortex-line reconnection (the analogue of fiux-line cut- 
ting), recently filmed by Paoletti et al,^^ who analyzed 
the trajectories of solid hydrogen tracers in the superfiuid 
to identify tens of thousands of individual reconnection 
events between quantized vortices. 

An extension of critical-state theory is needed to pro- 
vide the theoretical basis for calculating ac losses in 
superconducting power transmission cables fabricated 
from multiple helically wound layers of second-generation 
YBCO tapes. The helical currents generate longitudi- 
nal magnetic fields inside the cable, such that the super- 
current density J has components both perpendicular 
and parallel to B. 

In this paper, I extend critical-state theory to account 
for both flux transport and flux cutting in type-II su- 
perconductors. In Sec. II, I set down several of the basic 
equations needed and define the parallel and perpendicu- 
lar components of J and E. I use capital letters to denote 
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macroscopic fields of practical interest, J, and 

which in general depend on position r and time t. These 
fields are local averages over a length scale several times 
the characteristic mesoscopic lengths of type-II supercon- 
ductivity, namely the penetration depth A, the coh erence 
length ^, and the inter vortex spacing a ~ ^/(j)o/B. It has 
been known for over 50 years that in type-II supercon- 
ductors the spatial variation of the mesoscopic fields 6, 
jf, and e over these mesoscopic length scales is deter- 
mined by vortices, which carry magnetic fiux quantized in 
units of = /i/2e, the superconducting flux quantum.^ 
I will focus here on using an extended critical-state the- 
ory to calculate the magnetic-field and current-density 
distributions just above the critical current of a type-II 
superconducting cylinder, where J, and E are 

time-independent. However, since E is here produced by 
dynamic processes at the mesoscopic length scale involv- 
ing flux transport and flux cutting, we should keep in 
mind that 6, /i, j, and e are all time-dependent quanti- 
ties. 

In Sec. Ill, I review how to use the standard critical- 
state theory to examine B^ J, and E at the critical 
current Ic of the cylinder in zero applied magnetic field. 
In Sec. IV, I show how to use an extended critical-state 
theory to calculate the magnetic-field and current-density 
distributions at Ic in an applied longitudinal magnetic 
field Ha. In Sec. V, I calculate the dependence of Ic upon 
Ha when the critical current densities at the thresholds 
for fiux transport and fiux cutting depend upon the local 
value ofB. In Sec. VI, I show how an applied current den- 
sity affects the stability of a helical vortex arc stretched 
between two strong pinning centers. The results pro- 
vide us with a new model for the angular dependence of 
the critical current density, which simultaneously initi- 
ates fiux cutting and fiux transport. In Sec. VII, I discuss 
how the theoretical results qualitatively explain a variety 
of experiments, and I then turn to the issues of force-free 
configurations, the interactions between fiux cutting and 
fiux depinning, and the extensions needed to treat time- 
dependent problems. 



II. BASIC EQUATIONS 

To describe the behavior at the critical current of a 
long type-II cylindrical wire of radius R subjected to a 
parallel applied magnetic field i^^, let us assume that the 
local magnetic fiux density B = B{p)a{p)^ where a = 
sin a-\- z cos a, winds helically around the z axis with a 
pitch angle a{p), as shown in Fig. 1. Assume also that 
the local current density J winds helically around the z 
axis but at an angle (j){p) relative to B. The components 
of J along a and $ = a x p are Jy = J\\{p)a{p) and 
Jj_ = J±{p)${p)^ where /3 = 6 cos a — zs'ma. Similarly, 
when an electric field is generated, E = E\\a -\- E±$. 

In the steady state, when there is no time dependence 
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FIG. 1: Vectors used in this paper: The unit vectors in cylin- 
drical coordinates are p (out of the paper), (azimuthal), 
and z (longitudinal). The local magnetic induction B = Ba 
winds helically around the z axis with pitch angle a. The 
local current density J also winds helically around the z axis 
but at an angle (j) relative to B. The components of J along 
a and P — a x p are J\\ and J±. 

of B and J, Ampere's law and Faraday's law require that 
J = V xH, (1) 

V X = 0, (2) 

where H = {1/ po)\/ bF{B) and F{B) is the Helmholtz 
free energy density in the superconducting state. 



III. FLUX TRANSPORT AT THE CRITICAL 
CURRENT IN ZERO APPLIED LONGITUDINAL 
FIELD 

We assume that the superconducting wire contains a 
randomly distributed dense array of point pinning cen- 
ters. We also consider the case that the wire initially is 
in the fiux-free state, such that when a current / is ap- 
plied along the axis in the z direction in the absence of 
an applied magnetic field, azimuthal magnetic fiux pen- 
etrates in from the surface ai p = R with a distribution 
governed by a balance between the Lorentz force density 
J X B and the pinning force density Fp = pJc±B. In this 
case, wherever ^ ^ 0, a = 7r/2, d = 6>, /3 = — i, ^ = 0, 
J|| =0, and Jj_ = —Jj_z = Jc±z. We usually expect Jc± 
to be a monotonically decreasing function of 5, but for 
the present let us use the Bean model, for which Jc± is 
independent of B. 

If we carefully distinguish between B and p^H^ we 
note that since H{R) = //27ri?, no vortices can enter 
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the cylinder so long as / < /d, where Id = 27i RHd^ 
where Hd is the lower critical field. (We assume here 
that there are no surface barriers to vortex entry.) If 
Jci. < Hci/R^ the critical current is then /d, because 
any entering vortex ring of radius p will simply collapse to 
zero radius under its own line tension, which is (^oHd/p 
per unit length. 

When Jc_L > Hd/R, there is a minimum radius, 
Pc = Hd/Jc± < Ri at which a vortex ring can be held 
by the pinning forces. We seek an expression for the crit- 
ical current Ic at the first appearance of a steady-state 
longitudinal electric field, generated by the periodic nu- 
cleation of vortex rings at the surface p = R and their 
periodic self-annihilation when they are driven inwards to 
a radius slightly smaller than pc. According to Faraday's 
law, Eq. (2), £^ is a constant, independent of p, which is 
given by the Josephson relation E = B x v. Since here 
B = BqO and v = Vpp^ the electric field has only a z 
component, Ez = —BoVp. (Note that Vp = — since 
vortices are constantly moving inward.) 

When Id < I < Ic-, the solution of Eq. (1) is 
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Pci<P< R, (3) 



where 



Pel 



c± 



(4) 



is the radius at which H{p) = Hd and B{p) = 0. Note 
that Hci < H{p) < H{R) if pd < p < R. The critical 
current Ic is the value of / at which pd is reduced to pc, 
such that 



Ic = tt{R^ +pI)Jc±, 



Id 



2Rpc ' 



and 



H{p) 



ip' 



pI)Jc± 



2p 



pc< p< R. 



(5) 
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All the above results are greatly simplified if at the 
critical current pc = Hd/Jc± R^ which is equivalent 
to the condition that H{R) = Ic/2ttR ^ Hd- In high- 

superconductors it is well known that for H ^ Hd, 
B ^ pqH to good approximation. Thus if one makes 
the approximation that B = p^H from the outset and 
ignores line tension effects, this corresponds to setting 
Hd = and pc = in the above equations. When 
< / < /c, H{p) is then given by Eq. (3) but with pd = 
y/R'^ — I/7rJc± denoting the radius of the penetrating 
fiux front, where B = 0. At the critical current, which 
becomes simply Ic = 7rR^Jc±, the field distribution is 
then H{p) = pJc/'^, and the current density is Jz = Jc±- 
For simplicity, we shall assume in the rest of this paper 
that at the critical current the magnitude of H at the 
surface is much larger than Hd , such that the simplifying 
assumption B = pqH is a good approximation. 



IV. FLUX TRANSPORT AND FLUX CUTTING 
AT THE CRITICAL CURRENT IN AN APPLIED 
LONGITUDINAL FIELD 

A. Extending the elliptic critical-state model 

While the behavior at the critical current in zero ap- 
plied longitudinal magnetic field is described as above 
in terms of widely accepted critical state concepts, what 
happens at the critical current in a finite applied field is 
not yet well established theoretically. I present here a 
theoretical description that I believe describes the fun- 
damental physics of the behavior under these conditions. 

Let us assume that both an electrical current / and a 
magnetic field Ha are applied along the z direction, par- 
allel to the axis of a long type-II cylindrical wire of radius 
R. Under these conditions we can expect the magnetic 
induction B and the current density J to wind helically 
around the z axis as indicated in Fig. 1. We need a num- 
ber of equations to determine how the physical quanti- 
ties vary with the radial coordinate p. For simplicity, 
we make the simplifying approximation B = poH. In 
cylindrical coordinates, Eq. (1) yields the two equations 



J± 



/dH 



H sin a 



H 



\dp p 
/ da sin a cos 
\dp p 



(8) 
(9) 



For Ha = 0, we had Jy = and a = 7r/2, so that Eq. 
(9) was satisfied, and we had J± = —Jc±, such that the 
solution of Eq. (8) for H was given by Eq. (3). 

However, what determines the values of Jy and J± 
when Ha > 0? One model, which we proposed in Refs. 
21-25 and has been called the generalized double- critical- 
state model (GDCSM), assumed that both |J||| = J^y 
and |J_l| = Jc_l at the critical current. Here Jc± was 
identified as the magnitude of J± at the threshold of 
fiux transport (depinning), and Jc\\ was identified as the 
magnitude of Jy at the threshold of fiux cutting. Recent 
experiments,^^ however, have found that two of the pre- 
dictions of the GDCSM, a cusplike angular dependence of 
the critical-current density Jc and a sawtoothlike behav- 
ior of the direction of the electric field E just above the 
critical current, are not seen experimentally. The smooth 
angular dependence of Jc observed experimentally^^"^^ 
is in much better agreement with an elliptic critical- 
state model, introduced by Romero- Salazar and Perez- 
Rodriguez. However, the angular dependence of the 
electric field E for J just above Jc was found to require 
an extension of the elliptic critical-state model, to be de- 
scribed later. 

Figure 2 shows the ellipse representing the critical cur- 
rent density Jc(5, 0) in the J± — J\\ plane, according to 
the original elliptic critical-state model, where is 
the angle between J and B = Ba, and J± = Jsin^ and 
Jy = J cos (p. In the O zone inside the ellipse described 
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FIG. 2: Behavior of a vortex array as a function of the mag- 
nitude J and angle (j) of the current density J relative to 
the direction of the flux density B — Ba when flux-line cut- 
ting and depinning interact and the critical current density 
Jc{B,(l)) is given by the elliptic critical-state model of Refs. 
32-34. 
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or 
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Je(B,0) = l/. 



(10) 
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neither flux transport (depinning) nor flux cutting oc- 
curs {E± = and = 0). Flux transport, for which 
the vortices are depinned and \Ej_\ > 0, occurs every- 
where outside the ellipse (except when J± = 0) in zones 
with labels including the symbol {Ej_ > 0) or T_ 
{E± < 0). Flux-line cutting, for which |£^||| > 0, occurs 
everywhere outside the ellipse (except when Jy =0) in 
zones with labels including the symbol C+ {E\\ > 0) or 
C_ (£^|| < 0). In other words, except for the special cases 
of ^ = ±7r/2, 0, and tt, the critical current density for flux 
depinning and flux cutting is the same, Jc{B^(j)). When 
(j) = ±7r/2, flux transport (but no flux cutting) occurs 
when J > Jc^ and when = or tt, flux cutting (but no 
flux transport) occurs when J > Jc. 

The magnitude of the component of the critical current 
density along J_\_ associated with depinning (see Fig. 2) 
is 



For fixed Jy , the threshold for depinning is reached when 
\J_\_\ increases to Jed- The magnitude of the component 
of the critical current density along Jy associated with 
cutting is 



Jcc{B,(j)) = Jc{B , (j))\ COS ( 



(13) 



For fixed Jx, the threshold for cutting is reached when 
|J||| increases to Jcc- In contrast to the assumptions of 
the GDCSM, important new features of the underlying 
physics within the elliptic critical-state model are the as- 
sumptions that (a) the threshold Jed for depinning is re- 
duced as the magnitude of Jy increases, i.e., we can think 
of Jed is a monotonically decreasing function of the mag- 
nitude of Jy, and (b) the threshold Jcc for flux-line cut- 
ting is reduced as J± increases, i.e., we can think of Jcc 
is a monotonically decreasing function of the magnitude 
of Jx. 

I propose that the above relationships between the 
thresholds for depinning and flux cutting are an essen- 
tial general property of type-II superconductors, even if 
the mathematical form of the experimentally determined 
Jc{4>) deviates from the ellipse of Eq. (11) in a particu- 
lar material. In other words, (a) the threshold Jed for 
depinning is always reduced to zero as the magnitude of 
Jy increases to its maximum threshold value J^y and (b) 
the threshold Jcc for flux-line cutting is always reduced 
to zero as J± increases to its maximum threshold value 
Jc±. The underlying reason for this interrelationship is 
that when J is at an angle relative to the helical in- 
stability is the triggering mechanism that results in both 
depinning and flux cutting. A simple model calculation 
for the critical current for the helical instability of an iso- 
lated vortex subsequently leading to depinning and flux 
cutting is given in Sec. VI. 

In summary, at the critical current density of the el- 
liptic critical-state model, J = Jc{B,(j)), and J lies on 
the ellipse of Eq. (11), such that J±_ = Jc(5,^) sin ^, 
Jy = Jc{B^ (j)) cos 0, and 



Jc{B,4>) 



\J %\V? (j) - 



tan^ 



(14) 



. COS^ 



where tan(/)c = Jc±{B)/ Jc\\{B). Note that within this 
model Jc±{B) = Jc{B,7t/2) and Jc\\{B) = Jc(5,0). 

Additional relations are needed to connect the compo- 
nents of J and E. Here we extend the original elliptic 
critical-state theory^^^^ by introducing the general rela- 
tions 



E^ 
Ell 



P±J±, 
Pll^lh 



(15) 
(16) 



Jcd{B,(j)) = Jc(5,(/))| sin( 



(12) 



where and py are nonlinear effective resistivities, 
which are expected to depend upon the current density 
J and magnetic flux density B. Since the physics of flux 
cutting is different from that of flux depinning, we expect 
that p\\ is not the same as p^, and hence we do not ex- 
pect E to be parallel to J in general. As shown below, a 
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useful parameter is the ratio r = p\\/p±, which has been 
found experimentally to be independent of J just above 
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One could use the following model for the effective re- 
sistivities p± and p\\: 

^± = 0, < < Jed, (17) 

= Pd{J± - Jcd)^ J± > Jed, (18) 
^11 =0, < J|| < Jee, (19) 

= Pc{J\\ - Jcc), J\\ > Jcc, (20) 



such that E = ^ E]_ + and J = ^ J]_ + J| obey 

^ = 0, J<Jc{B,(t)), (21) 

= p[J- Je(5,0)], J>UB,c^), (22) 



where p = y p^sin (j) -\- p^ cos^ 0. Here, corresponds 
to the flux- flow resistivity of Kim, Hempstead, and 
Strnad^^ and of Bardeen and Stephen, and pc is a pos- 
tulated analogous quantity. (The subscripts d and c refer 
to depinning and cutting.) The model of Eqs. (19)-(22) is 
capable of describing the electric field components when 
there is a clear-cut linear onset of an electric field as J 
crosses the ellipse Jc(5, 0) shown in Fig. 2. On the other 
hand, if there are fiux-creep effects that make the resistive 
transition more rounded, causing the critical-current den- 
sity to depend upon an electric-field criterion, the more 
general model of Eqs. (15) and (16) should be used. 

We are now ready to write down all the equations 
needed to calculate the current- density and magnetic- 
field distributions at the critical current of a cylinder. 
Consistent with the usual critical-state approach for de- 
scribing the behavior in the absence of a longitudinal 
field, we assume that just at the critical current we 
have both \J±\ = Jed and |J||| = Jcc- For a current 
/ and a magnetic field Ha both applied in the z direc- 
tion we therefore have J± = Jc sin ^, where ^ < 0, and 
Jii = Jc cos 0, which results in the two equations 



Jc sin ( 
Jc cos ( 



/dH Hsin a 



\ dp 
/da 
\d^ 



^^/da since cos a 
H{ — + 



r 



(23) 
(24) 



An additional equation is needed to relate a and (/). 
This comes from Eq. (2) , which in cylindrical coordinates 
tells us that E = Eqz, where Eq > 0, so that E± = 
—Eq sin a and E\\ = Eq cos a. Combining these equations 
with Eqs. (15) and (16) yields 



tanc 



-r tana, 



(25) 



where r = p\\/p±. We see immediately for the special 
case of r = 1 that we would have (j) = —a, such that (see 
Fig. 1) J would have only a z component; Eq. (1) then 
tells us that the solution for H would have the property 
that Hz{p) = Ha- For arbitrary values of r, however, we 



can use Eq. (25) to eliminate (j) in favor of a, which yields 
the following first-order differential equations for H and 
a: 



dH _ H sin^ a 
dp 



rsma 



+ . 

P Jr'^ sin^ a/J^j_-\-cos'^ a/ 



da 
dp 



sm a cos a 



cos a 



P H^r'^ sin^ a/j2j_+cos2 



(26) 



,(27) 



I shah refer to the collection of Eqs. (10), (23), (24), 
and (25) as the extended elliptical critical- state models 
since the new equation, Eq. (25), goes beyond what was 
assumed in Refs. 32-34. Note, however, that these equa- 
tions are specialized for cylindrical geometry and would 
need generalization for other geometries. 



B. B consumption by flux cutting 

As was noted in Refs. 4, 8 and 37, the steady-state 
time-averaged voltage produced along a current- carrying 
type-II superconductor in a longitudinal magnetic field 
cannot be described as a fiux-fiow voltage generated by 
an inward collapsing array of helical vortices in the ab- 
sence of fiux cutting. The reason for this is that if 
no fiux cutting is occurring, the electric field is given 
hy E = B X t; ^5,35,38 B is the locally averaged 

magnetic fiux density generated by an array of vortices 
moving with a velocity v. According to Faraday's law, 
dB/dt = — V X a helical array of vortices continu- 
ously nucleating at the surface and moving inward with 
a velocity v while carrying a longitudinal component of 
B would produce an azimuthal component of the electric 
field, thereby leading to an ever-increasing longitudinal 
magnetic flux density. 

Flux cutting is the means by which the time derivative 
of B can be reduced to zero.^^ Multiplying Faraday's law 
by the unit vector a and making use of Eq. (9), we obtain 
the following equation describing the time rate of increase 
of the magnitude of in cylindrical geometry: 



dB _ dE± E^cos^a J\\E\\ 
~dt ~ dp p IT 



(28) 



The flrst two terms on the right-hand side of this equa- 
tion, which could be expressed as {dB /dt) transport, is 
simply the rate at which the local value of B is increased 
by the transport of B toward the cylinder axis with a 
velocity v = ExB / B'^ . For the experimental conditions 
considered here, {dB / dt) transport is always positive. The 
third term on the right-hand side, which could be ex- 
pressed as {dB / dt) cutting, is the rate at which the lo- 
cal value of B is increased as a consequence of flux-line 
cutting. However, since J\\E\\ is the rate of energy dis- 
sipation per unit volume, we see that {dB / dt) cutting is 
always negative wherever flux cutting is occurring. In 
the steady state, when dB/dt = 0, the local rate of in- 
crease of B due to flux transport is exactly balanced by 
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the local rate of decrease of B due to flux cutting, as can 
be shown with the help of Eqs. (23)- (25). 

According to the Josephson relation, in the presence of 
flux cutting, a uniform steady-state longitudinal electric 
fleld must be given by Eq = hv' /2e = 0o^^ where v' is 
the rate per unit length with which azimuthally directed 
flux quanta move inward and intersect a line parallel to 
the cylinder axis. Cyclic flux- cutting processes initiated 
by local helical instabilities^^^ allow helical vortices of 
one pitch to enter the cylinder and an equal number of 
vortices of longer pitch to exit the cylinder in such a 
way that on the average there is no net change in the 
number of longitudinal flux quanta in the cylinder. For 
this reason the time-averaged azimuthal component of 
the electric fleld is zero. 

We have seen in the above discussion that there is a 
competition between flux transport, which tends to in- 
crease B inside the cylinder, and flux cutting, which 
tends to decrease B. The efliciency of cutting relative 
to transport, as represented by the ratio r = p\\/p± [see 
Eqs. (15) and (16)] determines whether the average 
of Hz over the sample volume, is greater than, equal to, 
or less than the longitudinal applied fleld Ha. If r < 1, 
corresponding to a low efliciency of flux cutting, we ob- 
tain Hz > Ha and a paramagnetic longitudinal magnetic 
moment. If r = 1, corresponding to equal efliciencies 
of flux cutting and flux transport, we obtain Hz = Ha. 
If r > 1, corresponding to a high efficiency of flux cut- 
ting, we obtain Hz < Ha and a diamagnetic longitudinal 
magnetic moment. 
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FIG. 3: Solutions of Eqs. (26) and (27) vs p for r = 1, H{0) = 
1, and = 1. (a) i7 (solid), Hz (dashed), He (dot-dashed), 
and a (dotted), (b) Jz (solid), = (solid), Jy (dashed), 
and J± (dot-dashed) with the constants Jc\\ = 2 and Jc± = 1 
shown as dotted lines. 



C. Sample calculations 



1. No magnetic moment for r = 1 



I now present a few results of sample calculations of 
the magnetic-fleld and current- density distributions at 
the critical current of a type-II superconducting cylinder. 
Let Ha denote the longitudinal applied magnetic fleld, 
and assume that there is no surface barrier at the super- 
conductor's surface at p = i?, such that H{R) cos a{R) = 
Ha. (For simplicity, we consider here only the case 
Ha > 0.) By Ampere's law, at the critical current Ic 
where a flnite steady-state electric fleld E = Eqz flrst 
appears, we also have H{R) sma{R) = Ic/^ttR. The 
behavior depends in important ways upon the value of 
r = p\\/p± < 1 [see Eqs. (15) and (16)]. 

Depending upon the values of r (assumed here to be 
independent of J and H) and Ha^ the fleld along the axis 
H{0) = Hq can be either flnite or zero. When Hq > 0, 
we must have a{p) = k^p as p ^ 0, where = Jc||/2iJo, 
as can been seen from Eq. (27). On the other hand, if 
H{0) = 0, we may have a{0) = a^, where < < 7r/2, 
but then we must have H{p) = knp SiS p ^ 0, where 



k 



r sm ao 



(1 -h sin^ ceo) y^^^ sin^ <^o/-^cJ_ + ^^^^ (^o/Jc\\ 
as can be seen from Eq. (26). 



(29) 



Numerical solutions of Eqs. (26) and (27) for r = 1 are 
shown in Fig. 3. Note that Hz is a constant, equal to the 
applied longitudinal fleld when r = 1, as discussed in Sec. 
IV B. The reason for this is that when r = 1 [see Eqs. 
(15) and (16)], J is parallel to and since E = Eqz, J 
has only a z component, such that Jq = —dHz{p) /dp = 
and Hz = Ha. 



2. Paramagnetic moment for r < 1 

Solutions for r = 0.5 are shown in Fig. 4. When r < 1, 
J is no longer parallel to E^ and as fleld lines of H wind 
around the z axis as right-handed helices with a pitch 
angle a{p) > 0, streamlines of the current density J also 
wind around the z axis as right-handed helices with a 
pitch angle [see Eq. 25] 



■ a-\- ( 



tan 



(1 — r) tana 



1 + r tan^ a 



(30) 



Since tano^ > 0, t/^j > when r < 1. Since Jq = 
—dHz{p)/dp > 0, Hz decreases monotonically with in- 
creasing p, resulting in a paramagnetic magnetic moment 
per unit volume Mz = Hz — Ha > 0. As discussed in Sec. 
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FIG. 4: Solutions of Eqs. (26) and (27) vs p for r = 0.5, 

H{0) = 1, and = 1. (a) H (solid), (dashed). He 
(dot-dashed), and a (dotted), (b) Jz (solid), Je (solid), Jy 
(dashed), and J± (dot-dashed) with the constants Jc\\ = 2 
and Jc_L = 1 shown as dotted lines. 



IV B, B consumption due to flux cutting [see Eq. (28)] 
is needed to prevent an ever- increasing buildup of longi- 
tudinal flux. When r = p\\/ p± < 1, flux cutting is less 
efficient in consuming and this allows a larger value 
of Hz in the steady state. 

When r < 1 , the application of even a very small longi- 
tudinal magnetic field Ha can lead to a much larger value 
of the field Hz along the axis of the cylinder. Calculations 
illustrating this effect are shown in Fig. 5. 



3. Diamagnetic moment for r > 1 

Solutions for r = 2 are shown in Fig. 6. When r > 1, as 
field lines of H wind around the z axis as right-handed 
helices with a pitch angle a{p) > 0, streamlines of the 
current density J wind around the z axis as left-handed 
helices, i.e., with a negative pitch angle ipj < [see Eq. 
(30)]. Since Jq = —dHz{p)/dp < 0, Hz increases mono- 
tonically with increasing resulting in a diamagnetic 
magnetic moment per unit volume Mz = Hz— Ha < 0. As 
discussed in Sec. IV B, fiux cutting [see Eq. (28)] prevents 
an ever-increasing buildup of longitudinal flux. When 
r = p\\/p± > 1, flux cutting is more eflicient in consum- 
ing B and this results in a smaller value of Hz in the 
steady state. 




FIG. 5: Hz calculated from Eqs. (26) and (27) vs p for r = 0.5, 
Jc\\ =2, Jc_L = 1, and R = 1, showing a pronounced buildup 
of longitudinal flux along the cylinder axis even for very small 
values of the applied field Ha = Hz{R). 
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FIG. 6: Solutions of Eqs. (26) and (27) vs p for r = 2, i7(0) = 
1, and R = 1. (a) if (solid), Hz (dashed), He (dot-dashed), 
and a (dotted), (b) Jz (solid), Je (solid), Jy (dashed), and J±_ 
(dot-dashed) with the constants Jc\\ = 2 and Jc± = 1 shown 
as dotted lines. 



4- Critical current vs Ha 

Shown in Fig. 8 is the critical current ic = Ic/Ic±^ 
where Ic = 2ttRH0{R) and Ic = 7rR^Jc±, as a function 
of Ha = Hz{R), obtained Eqs. (26) and (27), vs p for 
several values of r assuming R = 1, Jc\\ = 2, Jc± = 1, 
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FIG. 7: Solutions of Eqs. (26) and (27) vs p for r = 2, J7 (0) = 
0, ao = 7r/4, and R = 1. (a) if (solid), ifj (dashed), He 
(dot-dashed), and a (dotted), (b) Jz (solid), Jg (solid), Jy 
(dashed), and J± (dot-dashed) with the constants Jc\\ = 2 
and Jcx = 1 shown as dotted lines. 
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FIG. 9: Jz [Eq. (31)] vs p for r = 0.3, 1, 3, and 10 with the 
constants J^w = 2 and Jc± = 1 shown as dotted lines. 



J|| cos a — Jl sin cc, and at the critical current we have 

r sin^ a + cos^ a 



J. 



^ r2 sin^ a/J^^+cos^ a/ 



(31) 



As shown in Fig. 9, Jz can be less than Jc± for some 
combinations of r and a, such that the integral over p 
yielding Ic can be less than /c_l, the effect being most 
pronounced for large values of r. 




FIG. 8: ic = Ic/Ic± (solid) vs Ha for r = 0.5, 1, and 2, 
calculated assuming = 1, Jc\\ — 2 (dotted), and Jc_l = 1 
(dotted), where /c_l = 7rR^Jc± (see text). 



and no 5 dependence of Jc (Bean model). 

Surprisingly, numerical calculations reveal that, de- 
pending upon the value of r, the critical current Ic = 
27r Jz{p)pdp can be less than Ici_ = 7rR'^Jc±. (For 
example, when r = 2, Ha = 0.113, ic = 0.9987, with 
all values of ic dipping slightly below 1 for small val- 
ues of Ha in Fig. 8.) The reason for this is that Jz = 



V. ACCOUNTING FOR FIELD DEPENDENCE 

OF Je_L AND Jell 

For simplicity, the above sample calculations were car- 
ried out using the assumption that Jc was a function of 
the angle (j) between J and B but was a constant indepen- 
dent of 5 or iJ, analogous to the commonly used Bean 
model. However, because the solutions of Eqs. (26) and 
(27) are obtained numerically, it is straightforward to in- 
corporate field dependence into Jc(5, ^). As an example, 
shown as the solid curves in Fig. 10 are sample plots of 
Ic for Ha applied parallel to the cylinder axis, calculated 
assuming that Jc{B^(j)) has the same dependence as in 
the Kim model, Jc{B) = Jc{0) / {1 -\- B / Bq) . For simplic- 
ity, in this paper we have assumed B = poH, so that for 
calculations of Ic when a longitudinal field is applied, J^y 
and Jc± in Eqs. (26) and (27) are replaced by their Kim- 
model analogs, Jc||/(1 + H/Hq) and Jc±/(1 + H/Hq). 
Note that self-field effects reduce the critical current in 
zero applied field to the value /c(0) = 2.47, below the 
Bean-model result, Ic± = t:E?Jc1_ = 3.14. Although Ic 
for parallel Ha depends upon r for relatively small fields, 
note that Ic ~ 7ri?^Jc||/(l + Ha/Ho) for large Ha, inde- 
pendent of r, because in this case a{p) and (/)(p) are both 
very close to zero for all p. 

When Ha is applied perpendicular to the cylinder axis, 
Ic{Ha) must reduce to its self- field value when Ha = 
and to Ic{Ha) ~ 7rR'^Jc±/{l + Ha/ Hq) when Ha is much 
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FIG. 10: Ic vs Ha parallel to / (solid) and Ic vs Ha per- 
pendicular to / (dashed) for r = 0.5, 1, and 2, calculated 
assuming R — 1, J^w — 10, Jc± — 1, and Kim- model field 
dependence of Jc with Hq = 1 (see text). The dotted curve 
shows 7ri^Vc||/(l + Ha/Ho), while 7vR^Jc±/{l + Ha/Ho) is 
indistinguishable from the dashed curve on this plot. 



greater than the self-field. However, application of an ar- 
bitrary perpendicular field destroys the azimuthal sym- 
metry of the field and current distributions, which now 
must depend on both the radial and azimuthal coordi- 
nates. Nevertheless, for simplicity, I have calculated the 
critical current using a = 7r/2 and H{p) = Hq{p) for 
p < R obtained as the solution of 



dH 
dp 



H 



P l + /ff^+flf/Ho' 



(32) 



an interpolation approximation that yields the correct Ic 
in both the self- field limit and the large- iJa limit. 



VI. SINGLE- VORTEX MODEL FOR Jc 

I next present a simple model that shows an inti- 
mate connection between the thresholds for depinning 
and flux cutting. Consider for simplicity a vortex seg- 
ment stretching between two very strong pinning centers 
at r± = (x^y^z) = (0, 0, =bc) (see Fig. 11). Assume that 
the vortex is subjected to a current density J = yJy-\-zJz, 
where Jy = Jsin^ and = Jcos^. (For simplicity, 
we here consider behavior in a Cartesian coordinate sys- 
tem rather than the cylindrical coordinate system used 
in Fig. 1. However, in both cases we may think of the z 
axis as the vortex direction, and (j) as the angle of J away 
from the vortex direction, as in Figs. 1 and 2.) Since the 
Lorentz force per unit length of magnitude J^^o due to 
the component Jy is perpendicular to the vortex line, it 
causes the vortex to bow out in the x direction, where 
the displaced vortex intersects the x axis at Xd- The 
bent vortex also experiences a Lorentz force due to the 
component Jz, causing the vortex to bend into the shape 
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FIG. 11: Helical vortex arc viewed looking down the x axis, 
calculated at the critical current for the helical instability 
when (j) = tt/S and c = 1 (the corresponding parameters 
are shown in the inset). The current density J is par- 
allel to the axis, and the helical arc wraps more than 
halfway around the cylinder of radius a (dashed), which is 
centered on the z' axis. The z' axis intersects the x axis at 
X = xo = Xdmax — ci = 0.201. The vortex arc intersects the x 
axis at Xdmax = 0.634, and the ends of the arc are attached 
to strong pins at r± = (0, 0, ±1). 



of a helical arc. For small values of J, the helical dis- 
tortion is stable, because the inward restoring force per 
unit length due to the vortex's line tension is able to 
balance the outward Lorentz force per unit length. How- 
ever, when J reaches a critical value, Jc, the helical arc 
becomes unstable and the vortex expands to ever larger 
radii, allowing it not only to escape the pinning centers 
and become depinned but also to grow outward, where it 
can meet other expanding vortices and undergo flux cut- 
ting. The critical current density Jc can be calculated as 
follows, using an extension of an approach used in Ref. 
39. 

Let us describe the helical arc using the coordinates 
{x\y\z') with unit vectors {x\y\z'), where the x' axis 
lies along the x axis and the y^ and z^ axes are tilted 
relative to the y and z axes by an angle (j) such that the 
z^ axis is parallel to J (see Fig. 11). Points on the hehcal 
arc are described by vector r' = x'a cos kz' — y'a sin kz' + 
z' z' . The origin of the primed frame is chosen to be at 
ro = (x^y^z) = (xo,0,0), where xq = Xd — a, such that 
the coordinates of any point in the primed frame are 
related to those in the unprimed frame via {x'^y'^z') = 
{x — xo^y cos <p — z sin 0, y sin (/) -\- z cos 0). In particular, 
the strong pins at r± = (0,0, =bc) have the coordinates 
r^^ = (xj_,2/±,^±) in the primed frame, and since these 
points lie at the ends of the helical arc, we have x^^ = 
a — Xd = acos(/cccos(/)), y'^ = =Fcsin(/) = sin^kc cos (p) , 
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and z'j^ = zbc cos ^. 

The unit vector tangent to the heUcal arc in the primed 
frame is 

, r' [z'^ —x'ka^\VLkz' — y'kaQ.o^kz'-\-z' 

1 [z ) — — 



ds' ViTW 
and its derivative is 
dV{z') _ N\z') _ k'^a{-x' cos kz' -\- sin kz') 



1 + Pa^ 



(33) 
, (34) 



where N'{z') = — x' cos kz' + y' sinkz' is the inward- 
directed unit normal to the hehcal arc and 



2^2 



1 + k^a 
k'^a 



(35) 



is the radius of curvature. 

The equation that determines the equihbrium radius a 
for a given current density is the force-balance equation 
fd^fr = where the inward restoring force is 



fr 



df'{z') eePaN' 



ds' 



1 + fc2a2 ' 



(36) 



and the outward driving force per unit length is, since 
J = Jz', 



fd = Jx cPof 



J(f>()ka 



N'. 



(37) 



The ratio of the magnitudes of the restoring and driving 
forces is 



R 



etk 



(38) 



and the two forces are exactly balanced when this ratio 
is equal to 1. 

In the absence of any pinning centers, Eq. (38) when 
a = tells us that a long, straight isolated vortex sub- 
jected to a parallel current density J is unstable to 
the growth of a helical perturbation of any wavevec- 
tor k < J^o/q, because then R < 1. For wavevec- 
tors k > J^o/q, on the other hand, it is possible for 
the restorin g and driving for ces to be balanced when 
a = au = ^(e^/J(/)o)^ — but this is a point of un- 

stable equilibrium. For a > a^^, < 1, the outward driv- 
ing force exceeds the inward restoring force, and the helix 
grows, while for a < a^^, > 1, the inward restoring force 
exceeds the outward driving force, and the helix shrinks 
to zero radius. Thus a straight vortex is stable against 
infinitesimal perturbations with wavevector k > J0o/e^. 
However, since fluctuations of all wavevectors k are possi- 
ble, we see that a long, straight isolated vortex is unstable 
to the growth of helical perturbations; in other words, in 
the absence of pinning, both J^w =0 and Jc± = 0. 

Returning to the model of the helical arc stretching 
between two strong pinning centers, we note that the pins 
place additional constraints on k and a via the equations 
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FIG. 12: / [solid, Eq. (41)] vs Xd = xa/c for (p = 0, 7r/32, 
7r/16, 7r/8, 7r/4, and 7r/2. For = tt/S, the forces are bal- 
anced {R = 1) for j = 1.4 (dashed) at two values of Xd/c 
(black points), a stable equilibrium point at Xds/c — 0.368 
and an unstable equilibrium point at Xdu/c = 1.080. As j 
increases, the two equilibrium points merge at the maximum 
of /, where j = jc and J = Jc. For = 7r/8, this occurs at 
Xdmax/c = 0.634 (large black point), where the dimensionless 
critical current density is jc{(t>) = 1.603. 



cos(A:ccos 0) = 1 — Xd/a and sin(A:ccos 0) = (c/a) sine 
Accordingly, we can write the force ratio as R{j, 0, Xd) • 
//j, where Xd = Xd/c and 



J 

f{(l>^Xd) 



J/ Jc±^ 

q/^oc, 

jka/^l^{kay, 



wsnr0-x|x 
\sm^6^x^J^^ 



cos 



ka 

= c/a = 2xd/{sm^ (/)-{- x^). 



COS( 



In the limits 



and 



(X), 



f{0,Xd) 
f{7r/2,Xd) 



Wl + (^^V2)^ 
2xd/{l^xl). 



(39) 
(40) 
(41) 

(42) 
(43) 



(44) 
(45) 



As illustrated in Fig. 12, when j is not too large, the 
restoring and driving forces are balanced {R = 1) at two 
displacements Xd of the vortex arc along the x axis, a 
stable point Xds and an unstable point Xdu- Note that 
for Xd slightly larger than Xds, R > 1 and the inward 
restoring force exceeds the outward driving force, such 
that Xd is driven back down to Xds^ but for Xd slightly 
smaller than Xds^ the opposite is true, R < 1 and Xd is 
driven back up to Xds- On the other hand, for Xd slightly 
larger than Xdu^ R < I and the outward driving force 
exceeds the inward restoring force, such that Xd is driven 
to ever larger values, while for Xd slightly smaller than 
Xdu, the opposite is true, R > 1 and Xd is driven down 
to Xds- As j increases to the value jc: the equilibrium 
points merge to the point Xdmax, where / has its maxi- 
mum value, /maa?(0) = jc ^ud R = 1. However, for any 
values of j > jc{(t^): the force ratio R = f /j is less than 



1, indicating that the outward driving force exceeds the 
inward restoring force for ah values of Xd, and the hehcal 
arc must expand to ever-increasing displacements. Thus 
when j > jc{(l>) or J > Jc{(fy^ the vortex arc undergoes a 
helical instability that not only leads to flux-line cutting 
with other vortices in the sample but also allows the free 
portion of the vortex to escape from the strong pins and 
undergo flux flow. 

Results for Jc and Xdmax vs (j) obtained using this sim- 
ple model are shown in Fig. 13. The key features are 
that (a) Jc|| is proportional to Jc_l, consistent with the 
experimentally observed correlation between these criti- 
cal currents, (b) the same helical instability leads to both 
depinning and flux cutting, showing that there is an in- 
timate connection between these processes, and (c) the 
threshold Jed for depinning is reduced to zero as the mag- 
nitude of J|| increases to its maximum threshold value Jc\\ 
while the threshold Jcc for flux-line cutting is reduced to 
zero as J± increases to its maximum threshold value Jc± . 

On the other hand, deficiencies of this model are that 
(a) the ratio Jc\\/Jci. = tt, but experimentally this ratio 
is apparently sample- and field-dependent, (b) it is an 
isolated-vortex model, which ignores the intervortex in- 
teractions that are generally important in type-II super- 
conductors, (c) it does not reproduce the experimentally 
observed angular dependence of JcC^), which is often 
well described using the elliptic model, and (d) the pin- 
ning model used here is greatly oversimplified and would 
need extensions to account for more realistic distribu- 
tions of pinning centers. The latter effect can be crudely 
estimated by convolving the helical-instability Jc(^) ob- 
tained above with a gaussian distribution function to ob- 
tain a pin-distribution average, 

Jc,avg{^) = [ g{(l)-^)JM)d^, (46) 

giij) = (l/^/^A</.)exp(-^VA<^'), (47) 

where is a measure of the width of the distribution 
of the pin-to-pin vectors around the z axis. Jc,avg{(t^) is 
plotted for Acj) = 7r/20 as the dot-dashed curve in Fig. 
13, where a corresponding plot of the elliptic model is 
shown as a dotted curve. 



VII. DISCUSSION 

A. Relation to experiments 

The theoretical calculations of Sec. JVC relate most 
closely to pioneering experiments by Walmsley.^ His ex- 
periments were carried out using a cylinder of Pb-40 at% 
Tl alloy, subjected to combinations of a magnetic field 
Ha and electrical current /, while both the axial mag- 
netization and the axial resistance were simultaneously 
monitored. As Ha was held fixed, a paramagnetic ax- 
ial magnetization was found to develop as the current / 
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FIG. 13: Jc(0)/Jc(7r/2) (solid) and Xdmax/c (dashed) vs 
(/)/(7r/2) obtained using the helical-instability calculation de- 
scribed in the text, for which Jc{7v/2) = Jc± = q/0oc 
and Jc(0) = Jc\\ = Tre^/^oc. The dot-dashed curve shows 
the pin-distribution average (46) obtained by convolving the 
helical-instability Jc(0) with the gaussian of Eq. (47) with 
A(f) = 7r/20, which yields Jc(0)/ Jc(7r/2) = 2.455. For compar- 
ison, the dotted curve shows the angular dependence of the 
elliptic model [Eq. (11)] for the same Jc\\/Jc± = 2.455. 



increased above the critical current Ic and flux flow oc- 
curred. This corresponds qualitatively to the behavior 
when r < 1 discussed in Sec. JVC 2 and shown in Fig. 4. 
The experimentally observed behavior can be interpreted 
as resulting from the nucleation of helical vortices pro- 
duced by the combination of the applied longitudinal field 
and the self-field. The helical vortices move toward the 
center of the cylinder, building up a longitudinal para- 
magnetic moment. As discussed in Sec. IV B, there is 
a balance between fiux transport and fiux cutting, but 
when r < 1, which corresponds to a low efficiency of ffux 
cutting for reduction of a paramagnetic longitudinal 
moment is produced. 

The experiments^ also revealed that when the applied 
magnetic field was reduced to a nominally zero value, a 
spontaneous longitudinal magnetic moment of either sign 
developed when I > Ic- Our interpretation of this effect 
is that, although the longitudinal field was adjusted to a 
very small value, it could not be made exactly zero, and 
as a result helical vortices penetrated from the surface 
during fiux fiow, carrying longitudinal magnetic flux to- 
ward the center. Figure 5 shows calculations for r < 1 
revealing large positive values of H^ in the interior of the 
sample even for very small positive values of Hz = Ha at 
the surface. Similarly, large negative values of Hz in the 
interior of the sample would also be expected theoreti- 
cally for very small negative values of Hz = Ha at the 
surface. 

The theoretical calculations of Sec. V provide a 
qualitative explanation of experiments by Karasik and 
Vereshchagin^^ showing striking differences between the 
critical currents of Ti-22 at% Nb and Ti-36 at% Nb wires 
measured in longitudinal and transverse magnetic fields. 
The critical currents in longitudinal fields Ha were all ob- 
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served to increase initially with Ha and to exhibit high 
peaks, where the longitudinal critical current was one or 
two orders of magnitude larger than the transverse crit- 
ical current at the same field. Such behavior can be un- 
derstood qualitatively as discussed in Sec. V and shown 
in Fig. 10. 

Various researchers experimentally studying the elec- 
tric field in longitudinal geometry have found inhomo- 
geneities in E along the sample length. ^'^^^^ However, 
such effects cannot be understood in terms of the present 
theory, which assumes uniformity of E just above the 
critical current. 



However, the above Bessel-function solutions are only 
one of infinitely many possible force-free solutions, and 
all such solutions fail to describe the full physics of the 
dynamical processes occurring at the critical current Ic 
of a superconducting cylinder in a parallel applied field. 
Force-free solutions do not satisfactorily explain the ori- 
gin of the longitudinal electric field that appears above 
Ic. If we had a truly force- free situation where J_l = 
and J|| > 0, then we would also have = and E\\ > 0, 
which violates Faraday's law. As a consequence, it would 
be impossible to balance the transport and cutting con- 
tributions to dB/dt appearing in Eq. (28). 



B. Force- free configurations 

In analyzing his experiments revealing a paramagnetic 
moment in current-carrying on Pb-40 at% Tl alloy cylin- 
ders subjected to a parallel magnetic field, Walmsley^ 
compared his results with force-free field theory, in which 
J is parallel to H. Using the terminology of the present 
paper, in cylindrical geometry force- free fields are those 
for which = 0. After setting = 0, we can integrate 
Eq. (8) to obtain"^^ 

H{p) = iJ(0)exp[- r r-^ sin^ a{r)dr], or (48) 

^0 

= H{R)exp[- / r-^ siB^ a{r)dr]. (49) 
Jp 

Since these equations hold for any arbitrary function 
a(p), there is an infinite number of possible force- free 
solutions. For each a{p), the corresponding J\\{p) can be 
calculated from Eq. (9). 

Walmsley^ focused his attention on Bessel-function so- 
lutions, which follow from the assumption that Jy (p) = 
kH{p), where [see Eq. (9)] 



da sin a cos a 
dp p 

is a constant. This equation is satisfied when^^ 
a{p) = tan"^[Ji(/cp)/Jo(/cp)], 



such that 



H,{p) = H{0)Jo{kp), 
He{p) = H{0)J,{kp), 

H{p) = H{^)^Jl{kp) + Jl{kp). 



(50) 

(51) 

(52) 
(53) 

(54) 



Since < H^{R) < H{0) and / = 27r RHo{R) > when 
< kR < 2.4048, solutions for k in this range yield 
paramagnetic force-free solutions for a cylinder of radius 
R carrying a current / in the z direction in a parallel field 
Ha = Hz{R). These solutions were found to be in good 
agreement with many of the measurements.^ 



C. Interactions between flux cutting and flux 
depinning 

Further theoretical work should be done to examine 
flux cutting in the presence of pinning centers. What 
initiates flux cutting evidently is the helical expansion 
instability first discussed in a geometries with linear di- 
mensions of the order of the penetration depth A,^'^^~-^^ 
such that interactions with the surface could stabilize the 
vortex against the helical expansion instability for small 
currents. However, Brandt^^'^^ showed that, because of 
their collective behavior, an array of vortices parallel to 
the axis of a cylinder of radius much larger than the in- 
tervortex spacing is unstable to a collective helical ex- 
pansion instability in arbitrarily small longitudinal cur- 
rents. These results suggest that pinning centers help 
to stabilize the vortex array against the helical expan- 
sion instability. This is consistent with the experimental 
observation that the longitudinal and transverse critical 
currents are roughly proportional to each other. 

The intimate relationship between flux-line cutting and 
flux depinning is also worthy of deeper study. Numerous 
experiments have revealed that the critical-current den- 
sities in longitudinal and transverse applied fields usually 
are roughly of the same order of magnitude, although the 
critical current in a longitudinal field is generally some- 
what higher than that in a transverse field. This suggests 
that the thresholds for flux-line cutting and flux depin- 
ning are closely linked. The model described in Sec. VI, 
in which the same helical instability simultaneously ini- 
tiates flux cutting and flux transport, should point the 
way to deeper understanding of this interrelationship. 



D. Accounting for both flux cutting and flux 
depinning in time-dependent problems 

While the present paper has dealt only with a steady- 
state problem in which there is no time dependence of 
the time- aver aged quantities and E, it should 

be possible to extend the above approach to solve prob- 
lems in which these quantities are time-dependent at fre- 
quencies of interest to power applications. The equations 
needed are simplest when the superconductor is macro- 
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scopically isotropic, i.e., when the penetration depth A 
and coherence length ^ are the same along different crys- 
tallographic directions, and the critical current densities 
for flux depinning and flux cutting depend only upon the 
magnitude of B and the angle (j) between J and B. The 
basic equations are then 

J = V xH, (55) 
V xE = -dB/dt, (56) 

where H = HH, H{B) = {1/ bF{B), and F{B) is 
the Helmholtz free energy density in the superconducting 
state. The displacement current can safely be neglected 
at low frequencies. H{B) can be obtained using standard 
magnetization measurements. Introducing the unit vec- 
tor a = B/B = H /H^ we can define the component of J 
along B as Jy = aJ\\. The component of J perpendicu- 
lar to B is then J± = J — J\\. Similarly, the components 
of E parallel and perpendicular to B are E\\= aE\\ and 
E±_ = E — E\\. What is also needed, but seldom experi- 
mentally determined to date, are the dependencies of p\\ 
and p_L in the expressions 

^11 = P\\J\\. (57) 
E^ = p^J^, (58) 

bearing mind that in type-II superconductors each of the 
quantities py and p± depends strongly upon its corre- 
sponding current density, with a well-defined increase 
when its threshold value {Jcc = Jc|cos(/)| or Jed = 
Jc I sin 01, as in Fig. 2) is exceeded. Finally, what is 
needed is a model for the behavior of Jc(B,0), similar 



to the elliptic critical- state model of Eq. (11), as shown 
in Fig. 2, where ([) is the angle between J and B. In many 
practical cases, it is likely that the time evolution of the 
magnetic-field and current-density profiles will need to 
be determined numerically by solving Eq. (56) step by 
step in time, as has been done in simpler cases in Refs. 
23 and 32. 

Further extensions of the above procedure would be 
necessary to incorporate the effects of anisotropy in both 
the intrinsic properties [i.e., if A, ^, and H{B) differ 
along different crystallographic directions] and extrinsic 
properties (e.g, if the pinning centers have linear charac- 
ter and are aligned along different crystallographic direc- 
tions) . 
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